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NATIONAL ADVISORY COMMITTEE FOR AERONAOTICS 
RESEARCH MEMORANDUM 

THEORETICAL ANALYSIS OF THE LONGITUDINAL BEHAVIOR OF 
AN AUTOMATICALLY CONTROLLED SUPERSONIC INTERCEPTOR 
DURING THE ATTACK PHASE 
By Ordway B. Gates, Jr., and C. H. Woodling 

SUMMARY 


A theoretical analysis has been made of the longitudinal behavior 
of an automatically controlled supersonic interceptor during the attack 
phase. The control system used to control the Interceptor's flight path 
was one in which a pitching velocity was commanded proportional to the 
longitudinal tracking error. Throughout the investigation the assumption 
is made that the target is flying on a straight-line path. 

Factors considered in this investigation included effects of control- 
system parameters, effects of limitations on control deflection and rate 
of control deflection, effects of initial tracking errors, effects of 
nonlinear variations in drag and lift with angle of attack and Mach num- 
ber, effects of nonlinear variations in pitching moment with angle of 
attack, effect of variations in interceptor forward velocity, and the 
effect of a nonnal acceleration limiter on the system performance. 

The control system considered in this investigation was .found to 
give acceptable control of the Interceptor ' s flight path dirrlng attack 
runs against a nonmaneuvering target. 

The inclusion of a nonlinear variation of drag and lift with angle 
of attack and Mach number resulted in relatively large variations in the 
Interceptor forward velocity during the attack rims. However, the effects 
of velocity changes on the overall responses during the attack runs were 
considerably reduced when a signal proportional to the change in forward 
velocity was fed back to the elevator servo. 


INTRODUCTION 


The Langley Laboratory of the National Advisory Committee For 
Aeronautics is presently engaged in an Interceptor research program, one 
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of the purposes of which Is'to evaluate the tracking performance of a 
supersonic interceptor equipped with various types of automatic control 
systems. The present paper is concerned with an analysis of that phase 
of the problem wherein the interceptor's radar locks on, with an initial 
vertical tracking error, to a bomber flying at a constant velocity; only 
maneuvers of the interceptor in the vertical plane are required to carry 
out the Interception. The results obtained from analysis of this longi- 
tudinal phase of the general tracking problem are Intended to provide 
information which will be useful in the synthesis of a satisfactory longi- 
tudinal control system for the interceptor being studied. The interceptor 
considered in this investigation is similar to that analyzed in reference 1, 
■which has a notched delta wing of aspect ratio 5*2 and 55° sweepback of 
the leading edge. 

For this investigation the Interceptor is assumed to be flying ini- 
tially in level flight at a Mach n\m±>er of 2.2 at an altitude of 
50,000 feet and the target is flying in level flight toward the inter- 
ceptor at a Mach number of 1.4, at various altitudes above 50,000 feet. 

No consideration was given to the effects of altitude changes on the 
interception problem discussed in this paper. 

The guidance equations presented in this paper are for a lead- 
collision type of navigation. 

The results of this investigation are presented, for the most part, 
in the form of interceptor and kinematic responses subsequent to radar 
lock-on, which were computed on the Reeves Electronic Analog Computer 
(REAC) . 


SmBOLS 

Iy moment of inertia about Y stability axis, slug-ft^ 

m mass of airplane, slugs 

c mean aerodynamic chord, ft 

S wing area, sq ft 

q dynamic pressure, Ib/sq ft 

V forward velocity, ft /sec 

M Mach number 


n 


nonmal acceleration, g units 
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g acceleration due to gravity, ft/sec^ 

q. 9c/2V when used as a subscript 

0 angle of pitch, radians unless otherwise specified 

a angle of attack, radians unless otherwise specified 

7 flight-path angle (7 = 0 - a,), radians unless otherwise specified 

u change in forward velocity, ft/sec 



relative change in forward velocity. 


u 

V 


elevator deflection, radians xmless otherwise specified 

-5e 


time, sec 



trim lift 

coefficient. 

Lift 

qS 

trim drag 

coefficient. 

Drag 


qS 


m 


pitching -moment coefficient. 


Pitching moment 
qSC 




7 } 


per radian 


per radian 




8Cm 


per radian 


per radian 


h 
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Cni 4 = ■: — —} per radian 
2Y 

SCjjj 

C|j)^ = , per radian 

2V 


Cm,,. = , per radian 


'D. 


d 


ha 


, per radian 


D differential operator^ — 

dt 

cf angle between interceptor X body axis and radar line of sight, 

positive when line of sight is above body axis, radians mless 
othenrwise specified 

R distance from interceptor to target along line of sight, measured 

positive from interceptor to target, ft 

a angular velocity of line of sight, (n = a + 0), radians/sec> 

positive when line of sight is rotating upward 

t(j time of flight of interceptor from instantaneous position to 

firing point, sec 

T time of flight of interceptor's rockets from firing point to 

predicted point of contact with target, sec 

M predicted miss distance, measured positive from interceptor to 

target, ft 

component of M along the instantaneoiis line of sight, positive 
when target is ahead of rockets at predicted time of impact, ft 

M|iu,q component of M perpendicular to the Instantaneous line of sight, 
positive when target is below rockets at predicted time of 
in5)act, ft 
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ey error in interceptor's flight path at any given instant, 

’’ Vjto + (Vj + Vb)t 

Tg elevator servo-system time constant, sec 

e^. output of filter, radians 

filter time constant, sec 

K tracking-loop gain constant, radlans/sec/radlan 

Kj. rate of pitch- feedback gain, radlans/sec/radlans/sec 

Kg elevator-servo gain constant, radians/radian/sec 

Subscripts ; 

I interceptor 

R rocket 

T target 

L limit 

ss steady state 

i input 

o initial value 


ANALYSIS 

Derivation of GixLdance Equations 


The type of navigation or Interception considered in this investi- 
gation is lead collision; that is, the Interceptor endeavors to fly a 
constant flight path such that at only one point on the path the rockets 
of the Interceptor may be fired and a hit obtained on the target. The 
rockets, subsequent to firing, fly a constant bearing course with the 
target to the predicted point of intact. The geometry of the attack 
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problem is shown diagrammatically in figure 1 . Generally, the vector 
eq.uation which must be satisfied is: 


E + — Vj-tQ + ^ t ^ + M ( l) 

The components of this vector equation along and normal to the instan- 
taneous line of sight are : 


E + Vy (t^ + cos + 0 - 7^^ = 
V^(^tQ + singer +9-7^^ = 




VjtG + ^Vj + Vj^^T 


cos (a + a) + Mje; 
sin((J + a) + 


> ( 2 ) 


The target flight-path angle 7,^ is taken as zero when the target is in 

level flight going away from the interceptor and is taken as jt when in 
level flight coming toward the interceptor. The equations may be rewritten 
in terms of the range, rate of change of range, and angular velocity of 
the line of sight as: 

E + E^tQ + - VpT cos(a + a) 

-EG^tg + - Vj^T sin(q + a) = 

where 

R = Vqi cos^ 0 + 0 - 7^j - Vj cos(0 + a) 

Rn = Vj sin(0 + a) - sin^0 + 0 - 7^^ 


= Mrs 


^NLS 


( 3 ) 
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The quantities a, Q, and a are defined as: 


R = Rq + / R dt 


a = Oq + / a dt = 0Q + / (a - e)dt 


e = 0o + / 0 dt 



In practice, R, a, and a would be available from the radar; however, 
in the analog solution of the problem these quantities were obtained from 
actual integrations. 

Certain simplifying assumptions were made in this investigation. 

The angles (cr + 9) and (c + a) were assumed to be small enough such that 
the cosines and sines of these angles are equal to unity and to the angle 
in radians, respectively. For these assumptions 


R 


cos 7,p + (cr + 0)sin 7,j, 


Rn « Vj(a + a) - 


{a + 0)cos 7^ - sin 7^ 


and the guidance equations are 




R + R^tQ + - VrT « M^s 

-Rfl ^t(j + rj - Vrt( 0 + a) « Mjjls 


> 


( 4 ) 
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The solution of equations (^i-) is accomplished hy computing continuously 
the value of (tQ + t) necessary for = 0 from the first of these 

equations and for this value of (t(j + t) computing the value of 

which will exist at (tg + t) seconds subsequent to the instantaneous 

time from the second of equations ( ^4-) . The time at which t(j is computed 

to be zero is the firing point for the interceptor's rockets. The time 
of flight of the interceptor's rockets is t, and throughout this investi- 
gation is assumed to be 1.5 seconds. The command to the control system 
is based on the error which exists at any time in the interceptor's 

flight path, which for this investigation was approximated by the 
expression 






+ (^1 + 


( 5 ) 


The validity of the foregoing assumptions in the guidance equations were 
checked by a digital solution on the Bell Telephone Laboratories X-66jkh 
relay computer at the Langley Laboratory where the exact guidance equa- 
tions were used. This comparison is discussed in a later section. 


Disciisslon of Flight -Path Control System 

The block diagram of the overall system is presented in figure 2. 
Briefly, the computed quantity is filtered, amplified, and used as 

the command to a pitch-rate command system. The dynamics of the filter 
and elevator servo are represented by simple first-order lag networks 

of the form i and i , respectively. The transfer function 

1 + 1 + rj) 

i is assimied to be representative of a low-pass filter which in 

1 + 

practice would be necessary to attenuate the high-frequency radar noise 
present in the computed command signal However, no attempt was made 

to include noise in the present investigation. The dynamics of the inter- 
ceptor were obtained from the linearized equations of longitudinal motion. 
For certain cases, these equations were modified to include specific non- 
linearities. All of the equations used in the analysis and the analog 
schematic diagram are presented in appendix A. The Interceptor parameters 
and other constants used in the analysis are presented in table I. The 
Interceptor stability derivatives and mass parameters were obtained from 
unpublished data and the results of reference 1. 
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SCOPE OF INVESTIGATION 


The Investigation may be conveniently divided into several parts, 
namely: 

(1) Preliminary determination of system gains. 

(2) Effects of nonlinear variations of drag, lift, and pitching 
moment on attack perfomnance. 

( 3 ) Effects of changes in forward velocity. 

(4) Study of the effect of system gains and other parameters. 

( 5 ) Effect of limits on deflection and rate of deflection of the 
elevator. 

(6) Discussion of normal acceleration limiter. 


EESULTS AND DISCUSSION 

Selection of System Gain Constants and Effects of Various 
Aerodynamic Parameters on Attack Performance 


Prior to the general study an investigation was made to determine 
values of the gain constants K and Kp for which the attack performance 

of the interceptor would be reasonably satisfactory as a starting point 
for the general study. The servo gain constant Kg is taken as unity 

throiaghout this paper. Examination of the block diagram (fig. 2) indi- 
cates that the forward-loop gain is KKg and the feedback gain is KpKg, 

and hence the assumption that Ks equals unity Imposes no restrictions 

on the system gain constants or perfomnance. From a theoretical analysis 
of the open-loop frequency response 7 /e^ and several preliminary rms 

on the REAC, acceptable values of K and Kp were found to be 5*0 
and 0.375> respectively. 

For these values of K and Kp tracking runs were coniputed for 

Ro = 60,000 feet, ^ = :tt, and do = 7 - 5 ° 45 ° smd are presented in 

figure 5* For these runs the interceptor motions were computed from the 
linearized equations of longitudinal motions presented in appendix A. 
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The responses shown in figure 5 include the predicted miss distance normal 
to the line of sight (Mijxs), interceptor normal acceleration (n), elevator 

deflection (5e)^ relative change in Interceptor forward velocity 



For these rims the elevator deflection was limited to ±20° 


and the rate of elevator deflection was limited to 'tl20°/sec. The miss 
distance My p-c» was started off-scale on the EEAC recorder in order to 
bring out more clearly the characteristics of in the vicinity of 

zero. For all runs presented in this paper, the transient responses are 
plotted up to the time at which the interceptor's rockets are assumed to 
be fired (tg = O). In view of the small angle assumptions made in deriving 

the guidance equations (eqs. (4)), the q-uantity (tg + t) Is dependent only 
upon the initial range Rq, target velocity rocket velocity Vr. 

rocket time of flight t, and Interceptor forward velocity Vj (l + u'); 

■*•0 

and for small values of u', the parameter (tg. + x) varies linearly with 
time. For Cq = 7-5° "th-e change in forward velocity is O.OTVIq, the 
value of MTjT.q a^t the assumed time of firing is -30 feet, and the peak 
normal acceleration is 7-6g. For Cq = 15° ^ the change in forward velocity 
is 0.l4Vj^, Mjjpg = -90 feet, and the peak acceleration is J.8g, For both 


values of Oq, the maximum perturbation in a was about 0.28 radian, but 
these transients are not presented. 


The miss distance Miyrr.q for neither of these runs is zero at the 

end of the run, but this condition is not due to the choice of gain con- 
stants. These nonzero values of Mi\rr.q can be attributed to the decrease 

in the Interceptor's forward velocity dirring the runs. The Interceptor 
is unable to maintain a condition of steady tracking ^n = 0, ~ 

as long as the forward velocity varies since, for the type of guidance 
considered, the interceptor must fly a constant flight path with constant 
velocity in order for the flight-path error e^ to be continuously zero. 

For these rims, u' never attains a constant value and consequently Mj^ 
is not zero at the assumed firing time. 


Effect of nonlinear variation of drag and lift with angle of attack 
and Mach number .- From unpublished wind-tunnel tests made for a model 
similar to the interceptor discussed in this paper, the variation of the 
drag coefficient Cp in the vicinity of the interceptor's trim angle of 

attack (oo = 0.033 radians) and initial Mach nimiber (Mq = 2.2) was found 



KACA RM L54K08 


11 


to be well approximated by the expression 


c„ = 0.014 

M V M 


il58„£ 

m 2 / 


and the variation of Cj^ with Mach mmber in this range was given by 


per radian 


If Mach number effects on CL and Ct are neglected, these expressions 
become, for Oq = 0.055 radians and Mq = 2 . 2 , 


Cp) = 0.027 + O.I 56 AL + 2.57(Z^)2 


Cp^ = 2.29 per radian 

The expressions for Cp) and Cp^, if only first-order changes in M 
considered, become 


(6) 


are 


Cp, = 0.027 + O.I56AX + 2.57 (Ai )2 


■^La 


2.29 I 



0.015 + 0 . 154 Zu + 2.05(Ax) 


2 ^ 
-I Mo 




> (T) 


Tracking responses were computed for K = 5*0, Kj. = 0.575^ 

Rq = 60,000 feet, and Gq = 7«5° 15° for the cases where equations (6) 

and (7) were used for Cp) and Cp, • The results are presented in fig- 

OG 

ures 4 (a) and 4 (b). Also shown in these figiires is the case for which 
ZiCp) is assumed to vary linearly with Ax, that is, AIIp, = O.I56AX. For 

Oq = 7.5° (fig. 4 (a)), the change in interceptor forward velocity is 
seen to be O.OTVj^ for the linear case, 0.15Vj^ when ADp, varies 
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nonlinear ly with Ax, and 0.l4Vj^ when AI!j) varies nonlinear ly with Ax 
and A1 and C]^ varies with Mach rnjmber. For the linear case 
Mjjig = -30 feet, and for the nonlinear cases -50 feet. As pointed 

out previously, the nonzero values of M]jls are due to the change in 
forward velocity. The same general trends are noted for cTq = 15° 

(fig. 4(b)), but the changes in forward velocity for this value of Oq 
are much greater for each drag condition investigated than were encountered 
for CTq = 7*5°* Also, the values of at the end of the run are 

larger for this value of Oq than for oTq = 7 . 5 °. 

Effect of feedback proportional to change in forward velocity .- A 
possible means of eliminating, or at least reducing, the value of M|^ 

at the assumed firing time is to feed back a signal to the elevator servo 
proportional to the change in forward velocity, such that a positive 
pitching moment is produced for a decrease in forward velocity. Rims 
are presented in figures 5(a) and 5(1^), for the case where Cj) varies 

nonlinear ly with Ax and A4, for Oq = 7*5° and 15 ° in which is incor- 
porated a feedback gain of 0.12. For each value of Oq the predicted 
value of Mjifxs is seen to be appreciably reduced. This type of feedback 
requires a bias error in the flight path in order for the interceptor to 
fly a constant flight path; but, on the basis of the predicted value 
of M]\[ls presented in figure 5 , this bias appears to be small. For com- 
parison runs in which the change in forward velocity was assumed to be 
zero are also shown in these figures. Despite the fact that the velocity 
changed considerably during these runs, there appears to be no appreciable 
difference between the cases which included the feedback proportional 
to u' in which forward velocity was allowed to vary, and the cases which 
neglected velocity changes. On the basis of these results, the remaining 
runs presented in this paper were computed with the assumption that veloc- 
ity changes can be made to have a negligible effect on the attack perform- 
ance of the interceptor being discussed. 

Effect of variations in the pitching moment due to angle of attack . - 
Recent results of wind-tunnel tests of complete models have often indi- 
cated a nonlinear variation of pitching moment with angle of attack. In 
order to check, at least qualitatively, the effect of a nonlinear pitching 
moment, runs were made for the assumed pitching variation presented in 
figure 6 (a) . This variation of Cm with angle of attack is generally 

similar to the type of variation obtained from wind-tunnel tests, but it 
should be pointed out that the range of a for which the pitching moment 
is nonlinear was arbitrarily selected on the basis of this range of a 
being in the range likely to be encountered in this particular problem. 

The results for this variation are presented in figure 6 (b) and afford 
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a comparison between the linear and nonlinear cases. For these runs 
K= 3 - 56 , Kp= 1 . 0 , 60,000 ft, and ffo = 7 * 5 °- altitude and 

Mach mamber of the Interceptor and target are the same as before. The 
most significant effects of the assumed nonlinearity appeared to be 
reflected In the Mi\rr.q » n, and 6 g transients and hence these are the 

variables presented In figure 6 (b). The maxlmim normal acceleration Is 
about 9 g for the nonlinear case as compared to approximately 7 S for the 
linear case. The nonlinear case also shows that there Is considerable 
overshoot In the transient, and the &g motion Is rather Irregular. 

It should be noted, however, that this assumed nonlinearity does not pre- 
vent the predicted from being approximately zero at the assumed 

time of firing. The results obtained for this assimied pitching -moment 
variation are in agreement with those which would have been Intuitively 
expected. The slope of the pitching-moment cvirve (fig. 6 (a)) Is seen 
to decrease in magnitude as the angle of attack increases , and finally 
reverses its sign. The general effect of reducing is to reduce 

the system spring constant in pitch and hence to increase the static 
sensitivity between pitching velocity (or normal acceleration) and elevator 
deflection. This effect Is reflected in the larger normal acceleration 
for the nonlinear case. 


General Effect of Control System Parameters and 

Initial Conditions on Attack Performance 

Tracking runs were computed which indicate the general effect on 
the system performance of variations in the gain constants and other 
system parameters . The comparisons are presented in the form of REAC 
runs and in the form of summary plots which indicate the effect of the 
various quantities on factors such as rise time, response time, and maxi- 
mum overshoot with respect to the miss distance Mfj-f.q , and maximum normal 

acceleration encountered dioring the run. In this investigation rise time 
is defined as the time required for to initially reach the point 

of zero miss and response time as the time for to reach and remain 

less than 30 feet. 

Effect of the tracking loop gain K .- The effect of the gain con- 
stant K on the system performance is shown in figure 7 values of 
0Q equal to 7 . 5 ° 15°. The results are presented for K = 

and 9- Foi’ these cases the initial range was 60,000 feet and the feed- 
back gain Kj. was equal to 0.375- The effect of K is seen to be 

essentially the same for both values of 0 q. As the gain is increased 
the system becomes more and more oscillatory and for K = 9-0 the Mjjjg 
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response is probably unsatisfactory for both values of Uq. The effects 

of K, with respect to rise time, response, time and overshoot of the 
Mjji 3 transient, and maximum normal acceleration encountered are summarized 

in figure 7(c) for Oq = 7*5° and 15°. As K is increased, there is only 
a slight variation in rise time and peak acceleration, which is due, to 
a large extent, to the fact that the elevator reaches its maximum deflec- 
tion of -20° very quickly and remains at that deflection for a time 
dependent upon K and Oq. As K Is increased, the response time tends 
first to decrease, and then become large as K is further increased. 

The response time for K = 9 and cTq = 15° is not shown, since for this 
combination of K and Oq, never reaches the condition where it 

remains less than 30 feet. For both values of Oq, the overshoot in 
Increases progressively with increases in K. 

Effect of rate feedback gain K^..- For K = 5*0, Oq = 7*5°7 and 
Eq = 80,000 feet, resxzlts are presented in figure 8(a) for values of 
K;^. = 0, 0.20, 0.375? 0.60, and 1.0. For Kj, = 0 the responses are seen 
to be rather oscillatory. As Kr is Increased the system becomes more 
stable, but the overshoot in is seen to increase with Increases 

in K^. These effects are due to the fact that, as the Interceptor is 

stabilized, its response to control Inputs becomes slower; hence the 
resultant increase in the Mi\rr.q overshoot. A summary of the effects 
of Kj. on the transient are presented in figure 8(b) . Inclusion 

of the pitch-rate feedback tends to reduce the response time and overshoot 
of and, in addition, to eliminate the oscillatory condition which 

exists for K^. = 0; but for values of greater than 0.375? "the over- 

shoot is larger than for Kj. = 0. There is seen to be a relatively neg- 
ligible effect of Kj. on the rise time which, as mentioned previously, 
is due to the fact that &g is at its limit of -20° at the beginning 

of the run and the rate feedback is ineffective during the early portion 
of the run. 

Effect of initial error in a .- For K = 3-0, K3. = 0-375? and 
Eq = 60,000 feet, several runs were made to evaluate the ability of the 
system to score a hit as the initial error in a is increased. The 
value of a actually reflects the altitude difference and the horizontal 
distance between the interceptor and target; at the beginning of the run 
the altitude difference is given by the expression Eq sln(ao + 0 q) . 

Eesults are presented in figure 9 for bg = 2°, 7-5°? 10°, 15°? and 20°. 

As Oq is increased, the time which elapses between reaching the correct 
flight path (Mjgxs = firing point becomes less and less, until 
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for 


Oo = 20° 


the maximum 

6e^ = -20°, 


the interceptor is unable to reduce to zero. The 

is a function of the range and 
In the present system for 


maximum allowable initial value of Or 


g the airplane can pull, 
nss is approximately 5 S- 


There is presented in figvire 10 a plot of the Initial range 
against the maximum initial cTq for which the present system could score 

a hit. This curve was obtained from a simultaneous solution of the 
equations 


Vl / cos 7 (t) 

° 0 L - 

= Pio cos(co + 0o) 


dt + ( Vj + Vj^ )t cos 


r(to) 


-V^ cos 7 ^(tQ + t) 


> (8) 






r(t) 

= Ro sin(cTo + 9o) 


dt + (Vi + Vj^jr sin 


7(to) 


- 1 ^ sin 7^(tQ + t) 


for the case where 7 ,^ = k. Equations ( 8 ) relate the horizontal and 

vertical distances traveled by the interceptor, rockets, and target to 
the horizontal and vertical distances which exist between the Interceptor 
and target at t = 0. If these equations were satisfied, a hit would be 
obtained. The function 7 (t) was calcixlated from the longitudinal equa- 
tions of motion for a step input on bg = - 20 °, which is taken as the maxl- 

mxjm value of 5e throughout this paper. Also presented in figure 10 
is the variation of I?q with cTq which was calculated on the assumption 
that 7 (t) = 7 ( 0 ) + The value of 7 used in this expression is the 

steady-state 7 due to an elevator deflection of - 20 °, and it was assumed 
that the airframe attained this output immediately upon application of 
control. The results obtained for this simplified approach to the problem 
are seen to be substantially the same as those obtained using the more 
exact approach. The curve indicates that for = 60,000 feet the 

largest Cq which can be used is approximately 22°. As a check, r\ms 
were made for a large value of K (in order to keep 5g = (Se)L for 
the entire run) and the results indicated that Oq = 20° was about the 
largest value which covild be tolerated. 
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In order to obtain a hit for Rq = 60,000 feet and cSq « 22°, the 

required time of flight of the interceptor from its initial position to 
the assimied firing point is seen from figure 10 to be l6.8 seconds. This 
explains the difference between the maximijm allowable Oq indicated by- 
figure 10 and the value obtained from the REAC results. For the KEAC 
rvins, the simplifications made in the guidance equations eliminated the 
dependence of (tg. + t) on Oq; and for the initial range Rq = 60,000 feet, 

t(j is equal to l4.8 seconds. Actually, 1 q for the REAC runs was 
slightly less than 1^4-. 8 seconds, since Rq was closer to 59,0,00 feet 

due to voltage limitations on the REAC. A check was also made for 
Rq = 30,000 feet and the result also agreed well with the curve of 
figiare 10. 

Effect of filter time constant Tp.- The effect of is Illustrated 

in figTjres 11(a) and 11(b) for = 0, 0.30, O.60, and 1.2 seconds. For 
these runs K = 3*0, Kj. = 0.375^ 7-5°^ and Rq = 60,000 feet. 

As is Increased from 0 to 1.2 seconds, the initial response in 

is seen to become progressively slower, which is due to the increased lag 
between the initial command €y and the response of the elevator motion. 

Also, there is seen to be an appreciable overshoot in the Mjjpg response 
for Tj. = 1.2 seconds, which is also due to the increased lag between 
and the elevator motion. The effects of T^. on the are 

svimmarized in figiire ll(c) for Oq = 2° and 7*5°* 

The values of Oq chosen to i l lustrate the effect of Tj have no 

special significance, but were used only on the basis that the results 
obtained for these values of Co were typical of the results obtained 
for all values of Oq up to the maximum allowable value for 
Ro = 60,000 feet. The same statement may be made concerning the values 
of Oq used in the subsequent sections . 

Effect of servo time constant Tg.- The effect of the servo time 
constant Tg is demonstrated in figure 12 for Tg = 0.03, 0.10, 0.20, 
and 0.30 second. The cases presented are for K = 3*0, Kj. = 0.375^ 

Oq = 2°, and Rq = 60,000 feet. As Tg is increased from 0.03 second 

to 0.30 second, the most important effect is that the overshoot in the 
Mfjis response is seen to become larger, but at the assumed firing point 

the miss distance is zero in either case. The effects of Tg on the 
Mjjls response are summarized in figure 12(b) . 
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Effect of control-surface limitations . - The effects of limitations 
on the rate of control-surface deflection and the magnitude of the sur- 
face deflection were investigated briefly for K = 3.0, Kj. = 0.375, 
and Rq = 60,000 feet. All of the results presented in this paper up 

to this point were obtained for the condition where 6 e was limited 
to ±120°/sec and 6 g to "i:20° . The effect of reducing the maximum value 
of 6 e to ± 6 o°/sec and 1 : 30 °/sec for this tracking run may be seen from 
flgiire 13 (a). For these runs Oq = 10°. The limiting control deflection 
was kept at ±20°. The effect on this run of reducing (Se)j^ was to reduce 
slightly the peak normal acceleration as (Se)^ reduced from + 120 °/sec 

to +30°/sec. This reduction in g causes a slight increase in the rise 
and response time of M| ^p , but it may be concluded that, at least for 

this case, the effects of limitations on the rate of control deflection 
were small. However, for cases where higher values of K would be 
required (for example, a maneuvering target) the effects of rate limita- 
tion would probably be much more important and shoiild be investigated 
thoroughly. The effects of limiting on the M^jj^ responses are 

summarized in figure 13 (t>) • 


Several runs were made for the case of (^e)]^ ~ il20°/sec and 
(&e) reduced from + 20 ° to ± 10 °; the results are presented in figure l4. 


For these cases Cq is equal to 5 °* The general effect of reducing the 
control deflection limits from 20 ° to 10 ° is to reduce the peak normal 
acceleration of the interceptor, and hence the rise and response time of 
the response. The maximum cTq for which M^c; can be reduced 

to zero is directly dependent on the limit imposed on 5g. 


Effect on tracking of limiting the command 9j^.- Since the result 

obtained from filtering and amplifying the error signal is used as 

the command to a pitch-rate command system, this command should be limited 
if it is desired to limit the Interceptor's normal acceleration from 
aerodynamic, pilot comfort, structiiral, or other considerations. The 
normal acceleration response to 0 ^ is given by 




n 


( 9 ) 
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Hence the steady-state normal acceleration can be limited to any desired 
value by limiting the input command 0^^. For the Interceptor being 

considered 



K 


s 



+ K^Ks 


1 

4.85 + Kr 


and, in general, 



Therefore, 


(6i)L 


(n)L(4.85 + Kr) 


_g_ 

Vl 


( 10 ) 


If 9^ is limited by this expression, it should be pointed out that 

only the steady-state n is being limited by use of this expression. 

The effectiveness of limiting the output transient acceleration^ by this 
method depends primarily on the response characteristics of (n/Si). If 

the system gains are chosen to give a fast response of n to 0^ with 
little or no overshoot, this means of limiting n should be satisfactory. 
The results presented in figure 15 afford a comparison of the cases for 
which there is no g-limiter and for the case where g is limited by 
equation (lO). For these cases K = 5 .O, Kp = 0.575^ Rq = 60,000 feet, 

CTq = 15 ° and (n)L = 5S* The airframe steady-state normal-acceleration 
response to the limiting value of 6g is approximately hence, the 

unlimited g case and the case of (n)j^ = 5 S should oscillate about the 
same value when &e = As can be seen, the peak n for the 

unlimited case is roughly 8 g. When (Qi)^ limited by equation (lO), 
the peak g response is reduced roughly to 6.5g. The rise time of MrjT.q 
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is seen to Increase slightly when 0^^ is limited, hut the response time 
is less than that for the imliinited case. 

Digital check on validity of simplified guidance equations . - In 
order to check the validity of the small-angle assumptions made in the 
guidance equations (see previous section entitled "Derivation of Guidance 
Equations'), a solution, using the exact guidance of kinematic equations, 
was obtained from the Bell Computer for comparison with the REAC tracking 
solutions which utilized the simplified equations. This comparison is 
presented in figure l6 for the case of Cq = 7*5°? Rq ~ 60,000 feet, 

K = 3-0, and Kj. = 0.375^ and the agreement is seen to be excellent. 


CONCLUSIONS 


The following conclusions were reached from a theoretical investi- 
gation of the longitudinal tracking behavior of an automatically con- 
trolled interceptor against a nonmaneuvering target: 

1. The control system considered in this Investigation (l.e., command 
on rate of pitch proportional to longitudinal tracking error) was found 

to give acceptable control of the interceptor's flight path during attack 
runs against a nonmaneuvering target. 

2. The inclusion of a nonlinear variation of drag- and lift with 
angle of attack and Mach number resulted in relatively large variations 
in the interceptor forward velocity during the attack runs. 

3- The changes in forward velocity computed for the runs in this 
investigation, although rather large, had a relatively small effect on 
the overall responses when a signal proportional to the change in forward 
velocity was fed back to the elevator servo. 

A. Consideration of a nonlinear variation of pitching moment with 
angle of attack which tended toward static instability at high angles 
of attack indicated that its primary effect was to increase the magnitudes 
of the interceptor's motions during the tracking runs. 

5 . The general effect of increasing the tracking gain K was to 
destabilize the tracking loop. 

6. Increases in the rate-feedback gain K^. tended to stabilize the 

Interceptor's longitudinal short-period oscillation, but had a destabili- 
zing effect on the tracking loop. 
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7 • The maximum initial angiiiarity between the interceptor ' s flight 
path and radar line of sight for which a hit can be obtained can be well 
approximated from the initial range and normal acceleration capabilities 
of the interceptor . 

8. Increases in either the filter time constant Tf or the servo 
time constant Ts had an adverse effect on the attack performance because 
of the increased lag between the input command and the elevator motions. 

9- The effects of limitations on the rate of control deflection did 
not appear to be large in this investigation for the limiting rates 
considered. 

10. The general effect of limitations on the magnitude of the con- 
trol deflection was to slow down the interceptor's responses, and hence 
the maximum initial tracking error which can be tolerated decreases as 
the limits on elevator deflection are reduced. 

11. Limiting of the interceptor normal acceleration was achieved 
by limiting the input command on pitching velocity, which effectively 
limits the interceptor's steady-state normal acceleration. 


Langlfey Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , October 22, 195^- 
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APPENDIX A 

EQUATIONS USED IN ANALOG OF TRACKING PROBLEM 
AND ANALOG SCHEMATIC DIAGRAMS 


All of the equations used in this investigation are presented in 
this appendix. The equations are presented in both symbolic and numerical 
form. It will be noted that certain of the equations presented in numer- 
ical form have been multiplied by constants. This was done in order to 
adjust the REAC voltages to satisfactory levels in the analog procedure. 

Linearized airframe equations (/j =0): 

■^o 


0 = Cn. — ^ — 9 + Cin- — ^ — a + C ^ a + 
2Vt. Iv ^ 2Vj^ ly “oc ly 


'lo -"Y 


p ' 3^ u' + C ' 8 ' 

Cmu ly ^ lY ® 


> (Al) 


i = e - Gj^ a - 2 Cl u' 


-'a mVi 


o 


mVr 


Clp, ' 6g ' 


HS . 


qS 


qS 


u' = -2 Ct 5 — u' - Ct — 0 - Ct, — a + Cr 


mVT hiVt 

•Lo -‘•o 


qS 

mVr 


a 
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It will be noted 
equations rather than 


that Cmgg * ) ^ 



and 6g ' appear in these 
. The relations exist that 



% 


This convention is adopted in order that the true physical phase relation- 
ships of the pitch-rate command loop be obtained for positive gain 
constants . 

Control Equations : 
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Kinematic eq^uations ( simplif iedj = it) 


to + T = 


Vrt - R 


R 


R = _v,p - Vi 


Rfi — Vj(n + <x) + V^(0 
MTjT.q = -Rfl(tG + t) - 

g 

Vi(tG + t) + Vj 


Vl = Vi^ + Vj^u' 


R = Rq + /R dt 


a = Oq + / a dt 


0 = 0n + / 0 it 


a 



(n - 0)dt 


a = 


m 


+ 0) 

'pT(a + a) 


2k 
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For the parameters presented in table I, the preceding eq^uations take 
the following form: 

Airframe eq^uations: 


0 = -0.18456 - 0. 01819a - 10.582a + 5 . 4696 e' 
= 0 _ 0.4565a - 0.0505U' + 0.05295e' 

-20ii' = 0.215u' + 0.5030 + 0.579a - 0.303a 


Control equations : 


2€f = 3.333e^ - 3-333ep 


01 = K^f 


0.055e' = 1.6670. - 1.667 Kj,9 - 1.667&„' 


Kinematic equations: 


tQ + T = 


3000 - R 

R 


12R = -16,320 - 12Vj 
20Rfl = 42,800(cr + a) + 27,200(a + 0) 


— 4Mjjx 2 = 4RJ2(tQ. + t) + 12,000(ct + a) 
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-20Vj = 
-R = 
-2a = 

-e = 
a = 

-loa = 


-42,800 - 42,800u’ 

-60,000 - f ^ d.'t 

-O.O66 - 2 a dt 

-0.035 - J'e dt 

Oq + ^ (n - 0)dt 

-lORn 

R 


The a and 0 appearing in equations (A3) are total a and 6, 
where'as the a and 0 In equations (Al) are perturbations away from 
the trimmed condition. The problem was slowed down such that 2 seconds 
of machine time was equivalent to 1 second of problem time. The scale 
factors used were: 100 volts = 100 seconds, 100 volts = 1 radian, and 

100 volts = 60,000 feet. 

In figures 17(a), 17(b), and 17(c) are presented the analog schematics 
of equations (Al), (A2) , and (A3)^ respectively. 
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TABLE I 

STABILITY DERIVATIVES AND MASS CHARACTERISTICS 
OF INTERCEPTOR AND OTHER CONSTANTS USED 
IN INVESTIGATION 


Altitude, ft . . 

p, slugs/cu ft . 
Vj^, ft /sec . . 

m, slugs . . . . 
ly, slugs -ft^ 

q, Ib/sq ft . . 

C , ft a • • • • 

S, sq ft a a a a 

Citiq^ per radian 
Cjn^, per radian 
Cm^, per radian 
Cjji^ ' > per radian 
Cl^, per radian 
C]3^, per radian 

cd 

Cl 

CmSg ’ radian 
CLg > per radian 

Vj^, ft/sec a a a 

V^, ft /sec a a a 


S0C « • • 

SGC • • • • 

Rq , f t a a a a 

0Q , radians 
oto ) radians 


. a a a 50,000 

a a 0 a 0003622 

2 , 11 + 0 (M = 2 a 2 ) 

a a a a 7 T 6 a 4 

a a 2 a 68 X 105 

826 

15 

Aoi 

a a a a - 2 a 8A 

a a a a -0 a 56 

a a a a ~0 a 2 8 

a a a a OaOO 

a a a a 2a29 

a a a a 0 a I56 

a a a a 0a027 

a a a a 0 a O76 

a a a a -Oa 295 

a a a a Oal65 

a a a a 2 , OOO 

a a a a 1 , 360 

a a a a 0a03 

a a a a 0 a 6O 

a a a a la 5 

a a a a 6O , 000 

a a a a 0,033 

a a a a 0a033 



Interceptor 

Figure 1.- Geometry of lead collision navigation used in present 
investigation for Mrg = 0. 
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Figure 2.- Block diagram of the longitudinal tracking system used in 

present investigation. 
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Figure 5 .- Interceptor and kinematic time histories for K = 5-0 and 

Kr = 0.575. 




Linear ACq , Linear C|_^ 

Nonlinear ACq CACq = f(Aa)3 , Linear C[_ 


Nonlinear ACoCACD=f(Aa,AM)D , 

Nonlinear Ci CCi =f(Aw)!) 
i-a '-a 

600 ■ 

. -20-^ 

« 300- 

^ -10 / V 


Figure ii-.- Effect of nonlinear drag and lift on interceptor attack 
performance. K = 5*0j K- = 0.375^ ~ 60,000 feet. 






g units ^^NLS > 


u' feedback excluded 
u' feedback included 


u' equal to zero 






(a) = 7 . 5 °. 


Figure 5 .- Effect of feedback, proportional to change in interceptor velocity 
on interceptor attack performance. K = 5 * 0 ; = 0 . 575 J I^o = 60,000 feet 



_2 I I I I I ■— 

0 2 4 6 8 10 

t, sec. 


J I I 1 I 1 1 t i _i I 1 ■ ■ i J 

12 14 16 18 0 2 4 6 8 10 12 14 16 18 

t, sec. 

(b) = 15°. 

Figure 5.- Concluded. 
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Pitching-moment coefficient , 



-8 


-4 


> I I 

0 4 8 

a , degrees 


I I I 

12 16 20 


(a) Assiomed variation of Cj^ with a. 

Figure 6.- Effect of nonlinear pitching-moment coefficient with angle of 
attack on interceptor attack performance. 
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o^, degrees or tg, seconds 

Figure 10.- Combinations of Rq, Oq, and for which a hit may be 

obtained in this investigation. = 9o = 
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Figure I5.- Con5>arison of u 
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(b) Control equations. 
Figure I7.- Continued. 
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(c) Kinematic equations. 
Figure I7.- Concluded. 









